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Abstract. We consider a model describing premixed combustion in the presence of fluid flow: 
reaction diffusion equation with passive advection and ignition type nonlinearity. What kinds of 
velocity profiles are capable of quenching (suppressing) any given flame, provided the velocity's 
amplitude is adequately large? Even for shear flows, the solution turns out to be surprisingly 
subtle. In this paper, we provide a sharp characterization of quenching for shear flows: the flow 
can quench any initial data if and only if the velocity profile does not have an interval larger 
than a certain critical size where it is identically constant. The efficiency of quenching depends 
strongly on the geometry and scaling of the flow. We discuss the cases of slowly and quickly 
varying flows, proving rigorously scaling laws that have been observed earlier in numerical 
experiments. The results require new estimates on the behavior of the solutions to advection- 
enhanced diffusion equation (also known as passive scalar in physical literature), a classical 
model describing a wealth of phenomena in nature. The technique involves probabilistic and 
PDE estimates, in particular applications of Malliavin calculus and central limit theorem for 
martingales. 

1. Introduction 

A mathematical model that describes a chemical reaction in a fluid is a system of two 
equations for concentration n and temperature T of the form 

T t + u-VT = kAT + Mg(T)n (1.1) 

n t + u ■ Vn = — An — Mg(T)n. 
Le 

Here k is the thermal diffusivity, Le the Lewis number (ratio of thermal and material diffusiv- 
ities), and M the reaction strength. The equations (jl.lj) are coupled to the reactive Navier- 
Stokes equations for the advection velocity u(x,y,t). Two assumptions are usually made to 
simplify the problem: the first is a constant density approximation 3J that allows to decouple 
the Navier-Stokes equations from the system (jl.lj) and to consider u(x,y,t) as a prescribed 
quantity that does not depend on T and n. The second assumption is that Le = 1 (equal 
thermal and material diffusivities). These two assumptions reduce the above system to a sin- 
gle scalar equation for the temperature T. We assume in addition that the advecting flow is 
unidirectional. Then the system (jl.lj) becomes 

T t + Au(y)T x = kAT + Mf(T) (1.2) 
T(0,x,y) = T (x,y) 
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with f(T) = g(T)(l — T). We are interested in strong advection, and accordingly have written 
the velocity as a product of the amplitude A and the profile u(y). In this paper we consider a 
nonlinearity / ^ of the ignition type 

(i) /(0) = /(l) = and f(T) is Lipschitz continuous on [0, 1], 

(ii) 36 > such that f(T) = for T G [0, O ], /(T) > for T G (0 O , 1), (1-3) 

(iii) f(T)<T. 

The last condition in (jl.Hj) is just a normalization. We consider the reaction-diffusion equation 
fll.2j) in the strip D = {i e I, y G [0, h}}. We take to be periodic with period h and with 
mean equal to zero: 

/ u(y)dy = 0. (1.4) 

A constant non-zero mean can be easily taken into account by translation. For the temperature, 
we impose periodic boundary conditions 

T(t,x,y) = T(t,x,y + h) (1.5) 

in y and decay in x. We will always assume that initial data Tq(x, y) is such that < Tq(x, y) < 
1. Then we have < T < 1 for alH > and (x, y) G D. For simplicity, we will usually assume 
that the initial data coincide with characteristic function of some set. More generally, we may 
assume that for some L and rj > we have 

T (x,y) > 6 + 7] for |x| < L/2, (1.6) 
T (x,y) = for \x\ > L. 

Equation (|1.2|) may be considered as a simple model of flame propagation in a fluid j2], 
advected by a shear (unidirectional) flow. The physical literature on the subject is vast, and 
we refer to the recent review ^H] for an extensive bibliography. The main physical effect of 
advection on front-like solutions is the speed-up of the flame propagation due to the large scale 
distortion of the front. The role of the advection term in (jl.2j) for the front-like initial data was 
also a subject of intensive mathematical scrutiny recently, see EI] for the references. 

The present paper considers advection effects for a different physically interesting situation, 
where initial data are compactly supported. In this case, two generic scenarios are possible. If 
the support of the initial data is large enough, then two fronts form and propagate in opposite 
directions. Fluid advection speeds up the propagation, accelerating the burning. However, if 
the support of the initial data is small, then the advection exposes the initial hot region to 
diffusion which cools it below the ignition temperature 8 Q , ultimately extinguishing the flame. 
The main purpose of this paper is to study the possibility of quenching of flames by strong 
fluid advection in the model (|1.2|) . The phenomena associated with flame quenching are of 
great interest for physical, astrophysical and engineering applications. For example, modelling 
of quenching and propagation of reaction fronts in fluid flow are relevant to studies of internal 
combustion engines, nuclear burning in stars and forest fires. Mathematically, the problem 
reduces to studying the advection-enhanced dissipation rate for the passive scalar. The passive 
scalar equation is one of the most studied PDE models, and has been the subject of extensive 
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research by both physicists and mathematicians. However, the question that we address here - 
controlling the rate of decay of the L°° norm in terms of the amplitude and geometric properties 
of the flow - while extremely natural, remained largely open until some very recent work. 

The problem of extinction and flame propagation in the mathematical model ()1.2|) was first 
studied by Kanel jHj in one dimension and with no advection. He showed that, in the absence 
of fluid motion, there exist two length scales Lq < L\ such that the flame becomes extinct for 
L < Lq, and propagates for L > L\. More precisely, he has shown that there exist Lq and L\ 
such that 

T(t, x, y) — > as t — > oo uniformly in D if L < Lq (1-7) 
T(t, x, y) — > 1 as t — > oo for all (x, y) 6 D if L > L\. 

We note in passing that it has been only very recently established by one of the authors that 
L = L\ in this situation [2111 • ^ n the absence of advection, the flame extinction is achieved 
by diffusion alone, given that the support of initial data is small compared to the scale of the 
laminar front width I = ■>/ k/M. However, in many applications quenching is the result of a 
strong wind, intense fluid motion, and operates on larger scales. Kanel's result was extended to 
non-zero advection by shear flows by Roquejoffre who has shown that (jl.7|) holds also for 
u^O with Lq and L\ depending, in particular, on A and u(y). However the interesting question 
about more explicit quantitative dependence of L Q , Li on A and u(y) remained completely open 
until recent work The following definition was given in 

Definition 1.1. We say that the profile u(y) is quenching if for any L and any initial data 
T (x, y) supported inside the interval [— L, L] x [0, h], there exists A = A (M, k, f,u, L) such 
that the solution of (|1.2j) becomes extinct: 

T(t, x, y) — > as t — > oo uniformly in D (1-8) 

for all \A\ > A . We call the profile u(y) strongly quenching if the critical amplitude of advection 
Aq satisfies Aq < CL for some constant C = C(M, k, f, u) (which has the dimension of inverse 
time) . 

The quenching property has been linked in jl] to hypoellipticity of a certain degenerate 
diffusion equation. In particular, one of the main results showed that u(y) is strongly quenching 
if there is no point y where all derivatives of u vanish. On the other hand, if u(y) has a plateau 
larger than a certain critical size, then u is not quenching. However hypoellipticity does not 
provide a precise solution of the problem at hand: a shear flow u{y) with a small plateau leads 
to an auxiliary equation which is not hypoelliptic, yet it is quenching. The first main result of 
this paper, Theorem 13.11 provides a sharp characterization of quenching shear flows. It states 
that a shear flow is quenching if and only if it has a plateau exceeding certain critical size. 
Nearby plateaux of smaller size will not lead to the same effect. This critical scale can be 
described in terms of existence of solutions to a nonlinear Dirichlet problem. The main new 
technical ingredient involves estimates on certain stochastic integrals, in particular application 
of Malliavin calculus to derive absolute continuity of the relevant random variables. 

The second goal is to study the dependence of quenching on scaling of the flow. Numerical 
experiments ^H] suggest that there is a certain scale of the flow for which quenching is most 
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efficient. Namely, if u(y) = sin ay, then the size La of initial data that can be quenched by flow 
Au(y) satisfies La ~ C a A with C a achieving maximum for some a^. Moreover, the constant 
C a satisfies C a ~ a -1 for large a and C a ~ a 2 for small a. Our Theorems 14.11 and 14.21 which 
apply to general shear flows, prove that in the small and large a asymptotic regimes one indeed 
has quenching for the initial data satisfying the above scaling. Central limit-type theorem for 
martingales is instrumental in obtaining the large a result. 

We mention that in a separate work jTH] one of us investigates the phenomenon of quenching 
in the presence of combustion- type reaction functions that do not have an ignition cutoff 8 , 
but are allowed to be positive for all T G (0,1). An important example of such function is 
the Arrhenius type reaction f(T) = e~ A ^ T (l — T), which is used in modeling of many chemical 
reactions. One of the main results of that paper, Theorem 1.3, is related to our Theorem 13.11 
It states that the quenching property of shear flows is in this case linked not only to the size of 
their plateaux but also to the decay rate of / at T = 0. Namely, if (in our setting) f(T) > cT p 
with p < 3, then no flow is quenching, whereas if f(T) < cT p with p > 3, then flows with small 
enough plateaux are quenching and those with large plateaux are not. This result, however, 
does not provide a sharp characterization of quenching flows, and should be understood as an 
extension to non-ignition reactions of the results of [3] rather than of our Theorem 13.11 

As a final remark we note that proving results for the system (jl.ljl is typically much harder 
than for a single equation ([1.2)1 . due to the lack of appropriate comparison principles. This is 
not so in our case. While in the paper we discuss quenching for a single equation ([1.2)1 . all our 
quenching results (including Theorems 13 . 1 f i) . 14 . 1 1 and !4.2j) extend immediately to the case of the 
system (jl.ljl . This is a consequence of a remnant of the maximum principle: the concentration 
n(t, x) remains bounded above by one for all times. Then T t + u ■ VT < kAT + Mg(T) in (Jl.ljl . 
and all the bounds we prove for quenching in the single equation model apply. 

The paper is organized as follows. In Section 2 we establish some auxiliary technical estimates 
on stochastic integrals. In Section 3 we prove results on quenching by shear flows and provide 
a characterization of the critical plateau size in terms of the corresponding Dirichlet problem. 
In Section 4 we deal with the scaling question. 



Results from this section will be used to obtain upper bounds on the solutions of (J1.2J) 
without the non-linear term, which can be expressed in terms of the Brownian motion. See the 
beginning of Section El for details and how this translates into estimates on the temperature T. 

We call a plateau of a function u G C(M) any maximal (w.r.t. inclusion) interval on which u 
is constant. We start by proving 

Lemma 2.1. Let u G C 1 (R) be bounded along with its first derivative and let Wjj denote the 
normalized one- dimensional Brownian motion starting at y. Then for any a el we have 
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Remarks. 1. In other words, the first probability is zero unless y is an interior point of a 
plateau of u with u{y) = |, in which case it equals the probability of W y staying inside this 
plateau for all s G [0, t]. 

2. This lemma for u G C°° and y not in a plateau of u follows from a probabilistic version of 
Hormander's theorem (see, e.g., Theorem 2.3.2]). Here we extend it to all u G C 1 and all 

y- 

3. We believe that the same result holds for u G C(K.) but we were unable to locate an 
appropriate reference in the literature. 

Proof. By Theorem 2.1.3 in JT] with F{W y ) = f Q u(W y )ds, the law of the random variable F 
is absolutely continuous with respect to the Lebesgue measure on M. whenever 



\DFf 2 = ) J q (J^ u'(W?)dr\ ds>0 (2.2) 



almost surely. We note that with the notation of |TTJ p. 24-26], if u G C 1 , then F G D 1,1 is the 
limit of F n (W y ) = i ELi w Wfc/n)' and DF ( S ) = f s u'(Wy)dr is the limit of 

fc=l fc=l j=fc 

Eq. ()2.2|) is obviously true if u' is not identically zero on an interval around y, that is, when 
y is not inside a plateau. In particular, for such y and all a, 

f(J^ u{W y s )ds = a^j = 0. (2.3) 

Now assume y to be inside a plateau J. For any open interval J with rational end points not 
intersecting any plateau of u, and any rational r G (0, t), let Bj T be the set of Brownian paths 
W y such that W y G J. Notice that every W v that exits / before time t, belongs to some such 

Bj :T - 

We have for any a 



P / u{W y s )ds 



W y G B J)T = 0. (2.4) 



This follows from ()2.3|) applied to the J* portion of the integral. Indeed — since U s = W. 



S + T 



(for s > 0) is just Brownian motion starting at W y , given any history {W / S S/ } S < T , the probability 
of J T u{W y )ds{= j Q ~ T u{U s )ds) being a - £ u(W y )ds is zero because W y is not in a plateau 
of u if W y G Bj^ T . By Fubini's theorem, ()2.4|) holds. Since there are only countably many sets 
Bj tT , the result follows. □ 

The main result of this section is 

Lemma 2.2. Let u G C 1 (1R) be periodic. Then for any compact interval S C (0, oo) we have 

P^jT u(W y )ds e [a,a + e] \ {tu(y)}^j ^0 (2.5) 
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as e — > 0, uniformly in (t,y,a) G S xE x R. 

Remarks. 1. Note that non-uniform convergence is an obvious consequence of Lemma 12. II 

2. The importance of this lemma lies in the fact that for large A it gives us a uniform (in 
(b,y,x) G S x R x R) estimate on the solution of lpOj) . (pO|) below, using (EHU1) . Through 
f)3.8p . ()3.9p this translates into an upper bound on the temperature T. 

To prove the lemma, consider the function 

p(t,y,a,s) = ^(WjOds G [a,o + e] \ {£w(y)} 

that is, the probability of u(Wj')£Zs G [a, a + e] and {u(W]?)} s <t not constant. 
Lemma 2.3. Under the conditions of Lemma \2. 6 A p is jointly continuous in R + x R 2 x 
Proof. For 6i G Rq and S 2 , 5 3 , 5 4 G R let 

6= ||«||oo|*i| + *IKI| oN + N + 



Then 

-i+<5i 

«(W7 + <y 2 )tfe- / u(W20<fe 







< Ikllool^ll + * 1 1 1 1 oo | ^2 1 



■ 



and we have 

|p(M-5i, y + S 2 , a + 8 3 , e + <5 4 ) - y, a, e) | 

<P^ u(W y s )ds G [a -5, a + 5] U [a + e - 5,a + e + 5} \ {tit(2/)}^ 

+P(exactly one of {w(Wj')} s < t and {^(WJ' + <5 2 )} s < t+( 5 1 is constant) . 
As 5 — > 0, the first probability goes to zero because by Lemma f2.1[ 

F (J u{W y s )ds G {a,a + e} \ {£u(y)}^ = 0. 

The second probability goes to zero because 

P({w(WJ)} s < 4 is constant) 

is continuous in (t,y). □ 

Proof of Lemma HOI By Lemma |2~H p(t,y,a,0) = 0. Hence by Lemma p(t,y,a,e) j as 
£ — > 0, for any (i, y, a). By joint continuity of p we then have p(t, y,a,e) J. as e — > 0, uniformly 
in (i, y, a) G -ft", for any compact -ft" C R + x R 2 . But p is periodic in y and p(t, y, a, e) = for 
M > £|M|oo + e. Thus p(t,y,a,e) j as e — > 0, uniformly in (t,y,a) G S 1 x R x R, for any 
compact S 1 C R + . □ 
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3. The Quenching Flows 

Let u(y) G C 1 (M) be a periodic function and let f(T) be an ignition- type non-linearity 
satisfying (i)-(iii) of ([1.3)1 . Let T(t,x,y), $(t,x,y), and ty(t, x, y) be the solutions of 

T t = K AT - Au{y)T x + Mf(T) (3.1) 

$t = kA$ - Au(y)$ z (3.2) 

*t = k* w - Au(j/)tf x (3.3) 

with (£, x,y) & ~Rq x M 2 and initial conditions 

T(0,:r,y) = $(0,:r,y) = #(0,x,y) = X[-l,l]( x )- ( 3 - 4 ) 

Notice that to prove quenching, one only needs to show 

||T(r, v )lloo<0o (3.5) 

for some r > 0. Indeed, the maximum principle then implies T(t, x, y) < 6 for all t > t. Hence 
we have 

T t = kAT - Au(y)T x (3.6) 

for t > t. To show (jl.Sj) we first notice that by integrating (j3.1j) in (x, y) G M x [0, h] (where h 
is the period of u) we have for || • || p = || • ||iP(R X [o,W) 

d 



9t -||T(t,.,.)||i = M / f(T)dxdy<M\\T(t, 

by ()1.3jl . and so ||T(r, •, •) ||i < oo. One can then, for instance, use the estimates on the parabolic 
kernel of the operator A — u ■ V for periodic divergence-free flow u on R n from JU] to show 
that the kernel for A — u ■ V on K x [0, h] with periodic boundary conditions is bounded above 
by Cr 1/2 for some C and all t > 1 (see [TH]). Therefore 

l|r(T + *, v )|| o<C7r 1 ^||r(T,.,.)lli (3-7) 

for t > 1 and (|1.8j) follows. Alternatively, there is a more elementary proof of ()3.7j) based on 
proving a Nash- type inequality for the evolution of (|3.6|) . namely that 

||r(r + t,-,-)|| 2 <(5t- 1/2 ||r(^-,-)lli 

with C7 independent of the flow. Such an estimate also leads to (|3.7|) by a duality argument. 
See [3 E] for more details. 

The functions $ can be used to estimate the non-linear evolution: 

T(t,x,y)<<S>(t,x,y)e Mt (3.8) 
sup $(t, x, y) < sup ^(t, x, y). (3.9) 

The first bound is achieved by replacing f(T) with T in (J3.1|) . while the second bound follows 
from the equality 

POO 

$(t,x,y)= / G(t,x - x')^{t,x',y)dx' 
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where 



G(x,t) = -^ =e - A /^t 



is the fundamental solution of the one-dimensional heat equation (using that \\G(t, = 1). 
The equality is verified by plugging it into (|3.2j) . 

Since $ and \& satisfy the above linear equations, we can apply the results from the previous 
section to obtain the following estimates. Let (W x ,W y ) be the normalized 2-dimensional 
Brownian motion starting at (x,y) and let (X x ,Y y ) be the random process starting at (x,y) 
and given by 

dX? = V2k~dW t x - Au{Y t y )dt, 
dY t y = y/2iidW?. 
Thus, Y t y = y + V2k~(W t y - y) = W y 2Kt and 

X x = x + V2~k{W? -x)- Au{Y y )ds = W» Kt - - / u(W*)da. 

Jo /K JO 

Then we have by (|3.2j) . ()3.4j) . and Lemma 7.8 in [TJ], 

/ A p2Kt 

<f>(t,x,y)=E(<f>(0,X?,Y?)) =W\W^- — J u{W y )ds G [-L, L] 
Similarly, 

*(t,x,y)=F\x-— u(W y )dse[-L,L}\. (3.10) 

The following result provides a sharp characterization of the quenching flows (see Defini- 
tion HH}. 

Theorem 3.1. With the above notation, there exists < £ < oo, depending only on M, k, and 
f , such that the following hold. 

(i) If the longest plateau of u is shorter than I, then u is quenching. 

(ii) If the longest plateau of u is longer than t, then u is not quenching. 
Moreover, this I is the infimum of all I such that the equation 

& = nA^ + Mf^) (3.11) 

on (x,y) G 1 x [0,/] with Dirichlet boundary conditions at y = 0,1, has a solution with 
0(0, •, •) compactly supported (and taking values in [0, 1]) such that does not go uniformly to 
zero as t — > oo. 

The key step in the proof is the following proposition. 

Proposition 3.2. For any l,L > let r(l, L) be the minimal time such that any solution of 
f|3.11|) on (x, y) e R X [0, 1] with Dirichlet boundary conditions at y = 0,1 and 0(0, •, •) supported 
in [—L,L] x [0,/] (and taking values in [0,1]), satisfies <fi(t,x,y) < 6q/2 fort > t(1,L). If such 
a time does not exist, we set r(l, L) = oo. Then with the above notation and I the length of the 
longest plateau of u we have the following. 
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(i) Ifr(l,L) < oo for every L < oo, then u is quenching. 

(ii) If t(1, Lq) = oo for some L Q < oo, then u is not quenching (and for all L > L and any A 
the temperature T(t,x,y) does not go uniformly to zero as t —> oo). 

Remark. Note that this result also applies in the case I = £. Whether quenching happens in 
this case depends not only on whether solutions of (j3.11j) with Dirichlet boundary conditions at 
y = 0, £, initially compactly supported, go uniformly to zero, but on this decay being uniform 
in all 0(0, •, •) supported in [— L, L] x [0,£] (for each L). 

Proof, (ii) Without loss of generality we can assume that the longest plateau of u is / = [0, 1]. 
Also without loss of generality, let u(0) = 0. Indeed — if u(0) ^ and T is the solution of 
()3.1|) with u(y) replaced by u(y) = u(y) — u(0), then T(t,x,y) = T(t,x — Au(0)t,y) and the 
result for T translates directly to T. 

Assume that for L = L the temperature T (with initial condition (|3.4|) and some A) goes 
uniformly to zero and let r < oo be such that ||T(r , ■, ■) ||oo < Oq/2. Since by comparison 
theorems any <ft from the statement of the proposition must satisfy <f)(x,y,t) < T(x,y,t) for 
all t, we obtain r(/, Lo) < To, a contradiction. Therefore T does not go uniformly to zero for 
L > Lq and any A, and so u is not quenching. 

(i) Fix L < oo, choose any < 5 < min{L, 1}, and let T2 = 1 + r(l, L + 5). We will show 
that for large enough \A\ and all x,y one has 

T(T 2 ,x,y)<6 . (3.12) 

This is (|3.5|) and so (|1.8|) will follow. 

Our strategy to show (|3.12j) will be to estimate T by $ via Lemma [2.21 for y outside of the 
plateaux of u (with large enough \A\), and by a suitable cf> from the statement of the proposition 
for y inside a plateau. 

Let d be the Lipschitz constant for /, define c = Mmaxjrf, 1}, and pick n G (0, 1) such that 



P( > £ ) < 7^- (313) 



2 J ~ 4e CT2 
This choice will become clear later. 

Let Aq be large enough so that for all \A\ > A , t 6 [n, t 2 ], x e M, and y not in the interior 
of a plateau of u we have 

(A I 2nt \ 9 

This is possible by Lemma f2. 21 and the first remark after Lemma f2.ll Then by (|3.8|) and (|3.9jl . 

for t G [ti, r 2 ] and y not in the interior of a plateau, 

supT(t,x,y) < e m sup*(t,x,y) < e Mt ^ < (3.15) 

In particular, f)3. 12j) holds when y is outside of all plateaux of u. 

We are left with the case of y inside a plateau. Hence consider a plateau J of maximal length 
I. All the following arguments will also apply to any other plateau of length I < I because 
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r(Z, L) < t(1, L) by comparison theorems. Therefore proving ()3.12|) for y G / will yield the 
same statement for all other plateaux of u and the proof will be finished. 

Again assume without loss of generality that / = [0,1] and u(0) =0. Increase A Q (if necessary) 
so that for any \A\ > Aq, 

"2kti 



Fl x- 



A 
2^ Jo 



u(W v s )ds G 



r 5 r 5 

L--.L+- 
2' 2 



< 



Oo 



4e CT2 



(3.16) 



whenever y G / and |x| > L + 5. Such A exists because by Lemma [2.21 



P 



2kti 9 k 

ui wv )ds G _ 



r 5 r 5 

L - - + x,L + - + x 



as 



|j4| — > oo, uniformly in y G / and x £ [-L — |,L+ |] (because w(0) = 0). Using (j3.8|) . 



(j3.13|) . and ()3.16|) . it follows that for y G / and \x\ > L + 5, 

T{r 1 ,x,y)<e MT ^{r 1 ,x,y) 



e Mn P( 



2,vn - 2^ 



4 /-2kti \ 

-/ e[-L,L] 



< e Mri P( x 



2k 



2kti 



w(W^)ds G 



L - -,L 
2' 



+ e 



9q_ 

4e CT2 



< 



0n 



(3.17) 

Next consider a function <p(t,x,y) defined on [ti,oo) x R x J, taking values in [0,1], and 
satisfying ()3.11)) with Dirichlet boundary conditions at y = 0, 1 and initial data 

#o 



T(n,x,y) 



< <f>(n,x,y) < X[-L-8,L+s](x)xi(y). 



2e c ( T 2- r i) 

Such a cf) exists because of f)3.15j) and f)3 . 1 7j) . and by the definition of r 2 we have 

(p(r 2 ,x,y) < y. 

Now let u = T - for (t, x, y) G [t x , oo) x R X /. Then 

w f = Aw + M[/(T) - /(#)] <Aw + c\oj\ 

By dSHHJ, 

w(ri,a:,y) < 

and by (l3~ToTl . 



^o 

2 e <-r 2 -T 1 ) 



sup{u(t, x, 0),u(t, x,l)} < 



0o 



2 e c ( T 2-ri) 

for t G [ti, r 2 ]. Now for Co = e~ ct u we have 

Cj t < Au; + c(|a)| — cD). 



(3.18) 

(3.19) 

(3.20) 
(3.21) 
(3.22) 
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Thus by ()3.21j) . ()3.22j) . and the maximum principle, 

0o 



2e CT2 

for t G [ti, r 2 ]. Thus u(r 2 , x, y) < 9 /2 whenever y G I. So by (|3.19j) . T(r 2 , x, y) < 9 for y G I. 
As mentioned before, this also holds for any other plateau of u. Together with (J3.15|) this gives 
()3.12|) . and the result follows. □ 

Proof of Theorem VJ . 1\ Let £ be defined as in the statement of Theorem 13.11 The fact that 
£ < oo is proved in jl] by constructing a subsolution of (|3.11|) on R x [0, 1] for large enough 
I. Proposition 13.31 below shows that £ > 0. Notice that by comparison theorems (see e.g. ^3 
Chapter 10]), a solution described in the statement of Theorem 13.11 exists when I > £ and 
does not exist when I < £. 

Let I be the length of the longest plateau of u and let r(/, L) be as in Proposition 13.21 
(ii) If I > £, then there exists a solution (f> described above. This means that \\(fi(t, -, > 6*0 
for all t. Since 0(0, x, y) is supported in [— L , -^o] x [0, 1] for some L < oo, we obtain r(Z, L ) = 
oo. Proposition I3.2n i) then gives (ii). 

(i) If I < £, then take 5 = (£ — l)/3. For any L < oo let be a solution of (J3.11j) on 
R x [—5, 1 + 5} with Dirichlet boundary conditions at y — —5, 1 + S and 

X[-L,L](x)X[0,l](y) < 0(0, X,y) < X[-L-l,L+l](x)X[-8,l+8](y)- 

Since I + 25 < £, there is r L < oo such that ||0(tl, -, ^Woo < ^o/2- By comparison theorems this 
implies t(1,L) < tl < oo. The result follows from Proposition !3.2lT j. □ 

In |3] an upper bound on i was provided by constructing a non-zero compactly supported 
<p(x, y) such that 

kA0 + M/(0) > 

in the sense of distributions. By comparison theorems, I is at most the diameter of the support 
of 0. Here we give a lower bound on £, in terms of the existence of a stationary ID solution of 

diiD). 

Proposition 3.3. Let I be the length of the shortest interval I such that there exists a non-zero 
if) : I —>■ [0,1], vanishing at the edges of I, such that inside I 

Kip" + Mf(ip) = 0. (3.23) 

Then £>£. 

Proof. Assume is a solution of (|3.11|) on (x, y) G R x [0, 1] with Dirichlet boundary conditions 
at y — 0, 1 and 0(0,-,-) compactly supported (and taking values in [0,1]), such that does 
not go uniformly to zero as t — > oo. Let be the solution of (|3.11|) with the same boundary 
conditions, but with 0(0, x, y) = sup x 0(0, x, y). By comparison theorems, > 0, and so also 
does not go uniformly to zero as t —>■ oo. 

Moreover, obviously 4>(t,xi,y) = 4>(t,X2,y) for any t,y,xi,x 2 , and so ip(t,y) = 4>(t,x,y) is 
well-defined and solves 

i> t = K^ n + Mf{i>). (3.24) 
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Since ip does not go uniformly to 0, Proposition 13.61 provides us if) solving ()3.23jh defined on 
[0, /]. A simple shooting argument can be used to prove that the set of all / for which solution of 
(I3.23|) does not exist is open. Thus the set of all I for which such solution exists has a minimum 
£, and the result follows. □ 

Corollary 3.4. With the above notation (and f(T) <T) we have £ > h^k/M. 

Remark. From results in [I] it follows that £ < c\Jk/M for some constant c depending on /. 
It follows that the critical plateau length £ is of the order of the laminar front width \J~k[M. 

Proof. By Proposition I3.3| there exists a solution if) of (|3.23|) on [0,£] vanishing at 0, 1. Since 
f(ip) < ip, we then have 

Kip" + Mip > 0. 

That is, the lowest eigenvalue of —A on [0,£] is at most Mk. Hence £ > ~KsJ k/M. But if 
£ = k^Jk/M, then necessarily ip{y) = csm(^jM/ Ky). This contradicts ()3.23|) because f{tp) = 
for small ■0. Hence £ > tt^k/M and Proposition 13.31 gives the result. □ 

The above provides also a criterion for strong quenching (recall Definition I l.lj) . 

Theorem 3.5. If the longest plateau of u is shorter than £ from Proposition \H. then u is 
strongly quenching. 

Proof. Going through the proof of Proposition 13.21 one observes that Aq depends on n, r%, and 
L. For L > 1 one can make T\ only depend on r 2 and from the conditions on A one then sees 
that as long as r 2 is bounded, A only depends on L. Moreover, this dependence is linear by 
Lemma 12.21 (and Lemma f2.1j) since the lengths of the intervals in f!3.14|) and (j3.16|) are O(L). 
Therefore as long as r(Z, L) is bounded in L, we have Aq < O(L) and so u is strongly quenching. 

Assume that the longest plateau of u is has length / and let 5 = (£ — 1)/3. Let ip solve ()3.24|) 
on [—5,1 + 6} with Dirichlet boundary conditions at y = —5,1 + 5 and ip(0,y) > X[o,i]{y)- By 
1 + 25 < £, the definition of £, and Proposition E21 ip goes uniformly to and so there is f < oo 
such that \\if)(f, -)||oo < #o/2- As in the proof of Proposition 13.31 comparison theorems show 
that if 4>(t, x, y) is any solution of (J3.ll)) on R x [0, /] with Dirichlet boundary conditions at 
y = 0,1 and taking values in [0,1], then <p(t,x,y) < ip(t,y). Hence t(1,L) < f for all L, and 
the result follows. □ 

The following proposition relates dynamical properties of reaction-diffusion equation with 
Dirichlet boundary conditions to existence of stationary solutions. Since we were not able to 
find this simple and natural result in the literature, we provide the proof in a slightly more 
general setting than needed for our application. Let Q be a bounded domain in M n with a 
smooth boundary. We also assume for the sake of simplicity that the reaction function is 
smooth. In one dimension this requirement can be removed and / only continuous is sufficient. 
This can be done by approximation from above with smooth /, comparison principles, and a 
simple ODE shooting argument. 
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Proposition 3.6. Assume that there is a solution of 

<p t = nA<p + Mf((p) (3.25) 

on (x,t) G Q x R + , with Dirichlet boundary conditions at dfl and </>(•, 0) compactly supported 
(and taking values in [0, 1}), such that <f> does not go uniformly to zero as t — > oo. Then there 
exists a positive solution if) : fl — > [0, 1] o/ 

kA^ + Mf{$) = (3.26) 

satisfying Dirichlet boundary conditions on dfl. 

Proof. For the sake of simplicity we let k — M — 1. Since by the maximum principle <f>(x, t) < 1 
for any t, standard regularity estimates imply that all Sobolev norms of (f>(x, t) are uniformly 
bounded in time: ||0(x, i)llff 3 (n) < C«- Define (f>-(x, t) = lim sup^^ (f>(x, t) at every x E fl. We 
claim that 0_ (x, t) is Lipshitz continuous and is moreover a weak subsolution, that is 

/ D(f)_(x)Dv(x) dx < / f(cj)_(x))v(x) dx 
Jn Jn 

for any v G Cg°(£)). To avoid certain degenerate cases, we define here Lipshitz continuity 
as \(f>_(x,t) — <f>-(y,t)\ < C\x — y\ for any x, y which belong to some ball B C fl, with the 
constant C independent of x, y and B. Indeed, let C\ be a uniform upper bound on \V(f)(x, t)\. 
Assume there exist x, y G B C fl with \(f>-(x) — 4>-{y)\ > 2C\\x — y\. From the definition 
of (/>_ it follows that there exist t n — > oo such that either <f)(y,t n ) — 4>-{x) > 2C\\x — y\ or 
4>(x, t n ) — <f)-{y) > 2Ci\x — y\. But this implies that for any e > 0, for all sufficiently large n we 
have \(f>(y,t n ) — (f>(x,t n )\ > 2C±\x — y\ — e, which contradicts the bound on the gradient of 0. 

Notice also that compactness of fl and uniform boundedness of |V0| show that 0_ is not 
identically zero and vanishes on dfl. 

Define A$(f)-(x) = 5~ 2 ^™ = i{4>-{x + $ e j) + <fi-{x — Sej) — 20_(x)), where e,- are unit vectors 
in coordinate directions. Next, we claim that for any x such that dist(x, dfl) > 5, we have 
As<p-(x) > —f((f)-(x)) — 7(5), where 7(5) converges to zero when 5 goes to zero. Indeed, by 
definition of (f>-(x), we have that for any e > 0, there exists a sequence t n — > 00 such that 
\<f>-(x) — <f>(x,t n )\ < e and (f>-(y) > (f)(y,t n ) — e for any y. Moreover, we can choose t n so that 
\<f> t (x,t n )\ < e. Now 

A S (j)_(x) = 5~ 2 J^(0_(a; + Sej) + (j)-(x - Sej) - 2(f)_(x)) 
3 

> -Ce5~ 2 + (T 2 ^2(<p(x + Sej, t n ) + <f>(x - be h t n ) - 2<p(x, t n )). 
3 

Using the mean value theorem and uniform upper bounds on derivatives of 0, it is not hard to 
show that 

(T 2 ^(0(x + 5ej,t n ) +(f)(x-5e j ,t n ) -2<f>(x,t n )) (fr XjXj (x,t n ) 
3 

uniformly in x and t n as 5 — > 0, with an error bounded by C5. Therefore, 

A 5 0_(x) > -Ce5~ 2 -C5 + Acf)(x, t n ) 
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> -f{<P{x,t n ))-C{e5- 2 + 5)-t 
>-f^_( x ))-C(e5- 2 + 5 + e). 

Since e is arbitrary, this leads to A$(p^(x) > —/(</>_ (a;)) — C5. 

Given v G C^°(f2), v > 0, such that dist(supp(w ), dQ) > 5, we have 

— / Ag<f)-(x)v(x) dx < / f((j)-(x))v(x)dx + C5\\v\\ L i(p i y 
Jn Jn 

Carrying out discrete integration by parts on the left hand side and passing to the limit 5^0, 
we get 

/ D<p_(x)Dv(x) dx < / f((p-(x))v(x) dx. 
Jn Jn 

Passage to the limit is justified since we know that <p-{x) is Lipshitz and therefore belongs to 

the Sobolev space W 1 ' 00 . Thus we see that 4>-{x) is a weak subsolution of ()3.2fi|l . 

Now consider initial data <p(x, 0) such that (j)-{x) < <p(x, 0) < 1. By the maximum principle, 

for all t we have <f)(x,t) > <p-(x). Consider 4>+(x) = liminf^oo (j>(x, t) > <ft-(x). By repeating 

the same arguments as above, we find that <p+(x) is a weak super solution. Then by well-known 

results (see e.g. jH], Theorem 9.3.1), there exists a weak solution if)(x) of (JS2EJ), satisfying 

4>-(x) < ip(x) < (j) + (x). By boundary regularity results, ip(x) is regular on all of Q. □ 

Results in this section extend without change to the case of shear flows in higher dimensions. 
The proofs are identical to those above, this time using higher dimensional Brownian motion. 
Assume that T(t,x,y) is a solution of (j3~Tj) . (EOJ) on l| x 1 x 1" with the C 1 shear flow 
u satisfying u(y) = u(y + hj6j) for j = 1,2, ... ,n and some hj > ({e 1; . . . , e n } being the 
standard basis in IR n ). 

The definition of quenching flows is identical to that for n — 1. A plateau of u is any maximal 
domain Q on which u is constant. We also say that a domain Q C R n is quenching if for every 
L < oo there is t(Q, L) < oo such that any solution <fi of (J3.11|) on (x, y) G IR x Q with Dirichlet 
boundary conditions on dQ and 0(0, ■, ■) supported in [-L, L] x Q (and taking values in [0, 1]), 
satisfies <p(t,x,y) < 9 /2 for t > r(0, L). Of course, the quenching property again depends on 
M, k, f . We then have 

Theorem 3.7. With the above notation the following hold. 

(i) // every plateau of u is quenching, then u is quenching. 

(ii) If u has a plateau that is not quenching, then u is not quenching. 

(iii) If dQ and f are smooth, and there is no non-zero ip : O — > [0, 1] satisfying 

nAtp + Mf(ip) = 

and vanishing on dQ, then the domain Q is quenching. Moreover, if each plateau of u is 
contained in some such domain, then u is strongly quenching. 

Remark. Note that if n > 2, then even non-constant u can have unbounded plateaux. 

Finally we note that we only considered initial conditions ()3.4j) for the sake of simplicity of 
presentation. It is obvious that our results apply also in the case of smooth initial conditions 
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satisfying, for instance, 

X[-l,l](x) < T(0,x,y) < X{~L-c L ,L+c L ](x). 

If we wish to consider initial temperatures that are not maximal (but still above the ignition 
temperature 9q) on an increasing family of regions, for example, 

VX[-L,L](x) < T(0,X,y) < 7]X[-L-c L ,L+c L ]{x). 

for some rj G (6q,1), then there is only one change — I in Theorem 13.11 is defined in terms 
of Dirichlet solutions <fi initially compactly supported and initially bounded above by rj. The 
above method actually applies in the case of any family of compactly supported initial conditions 
T L (0, x, y) as long as these are such that for any L\ and Si > there are L 2 and S 2 > so that 
T L2 {0, x 2 ,y 2 ) > T Ll (0,x 1 ,yi) - Si whenever \{x 2 ,y 2 ) - (an, 3/1) | < 5 2 (in particular, T L (0,-,-) 
continuous will do). This last condition is necessary for our proof of part (i) of Theorem 13.11 
because now we have 

-2nt 



$(t, x,y)=E (t l (o, < W s) ds i 



2Kt 

Here £ is defined in terms of (ft initially bounded above by the T L 's. 

4. Scaling 

In this section we study the dependence of the "quenching amplitude" , that is, the infimum 
of all A such that initial temperature distribution 

T(p,x,y) = xi-L,L](x) (4.1) 

leads to quenching, on the scaling of the profile of the shear flow u. Hence we consider 

T t = kAT - Au{ay)T x + Mf(T) (4.2) 

with u periodic and a > 0. The results of this section are motivated by and agree well with 
numerical simulations performed in ^Hl- The first is 

Theorem 4.1. Let u G C(R) be a periodic function with period h. Then there is C > 
such that for large enough a and \A\ > CaL, the solution of (|4.2j) with initial condition (|4.1|) 
satisfies T(t,x,y) — > as t — *> 00, uniformly in M 2 . 

Remark. The necessity of this bound can be explained by the fact that fast oscillations in the 
advection homogenize propagation of the flame (w.r.t. y) and so larger advection amplitudes 
are needed to expose the hot region to diffusion. 

Proof. Notice that we have 

sup^f— ,x,y \ = suprfx - — I u(aWg)ds G [-L,L]\ 
x \2k J x \ 2k Jo J 

(r 1 r 4kL 

/ u(aW%)ds G a, a + — — 
Jo L 1^1 
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supP( — / u(W? y )ds G 



supP 



u(W* y )ds e 



a, a + 



a, a + 



AkL 
W. 

LAI 



Let us estimate the last integral. 

First, we can assume u(y)dy = 0, since, as before, changing u by a constant does not 
change the result. Second, let v(y) be such that v'(y) = u(y) and v(y)dy = 0, and define 
z(y) = Jl v(s)ds. Hence, all three functions are periodic with period h. 

Now by the Ito formula (see, e.g., [TTJ Proposition 1.1.4]), 



v(W y )dW y + - / u(W y )ds 



almost surely. Thus, 



1 

a 



z(W y ) - z(y) 



u ( W y)ds = -(z(W« 2 ) - z{y)) - 2M(y, a, W y ) 



with 



Therefore with c 



M(y,a,W y ) = 



a 



v{W y s )dW y s . 



we have 



sup#( — ,x,y ) < supP( M(ay,a, W ay ) G 



2naL 4c 
a,a-\ — — - — I 



LA 



a 



From (|3.8|) and (|3.9j) we can see that to obtain (|3.5j) for r = (2k) 1 (and hence ()1.8|l ). we only 
need to prove 

2k Ac 

a, a + — H 

C a 



for some C and all large enough a. That is, 



< 



-Af/2/s 



swp¥(M(y,a,W y ) e [a,a + e]) < 9 e 



-M/2k 



(4.3) 



il-a 



for small e and large a. However, for each y, the family aM.(y, a, W y ) is a martingale with 
respect to a. It is not difficult to check that the central limit theorem for martingales (see, e.g. 
[HI, Theorem 7.7.3, or |14j ) applies to Ai(y,a,W y ) giving convergence in distribution to the 
normal random variable with variance 



a' 



h 



E 



v(w z s yds 



dz 



h 



\v(z)\ 2 dz > 0, 



where E denotes expectation with respect to the Brownian motion starting at z. Moreover the 
convergence can be shown to be uniform in y since all the estimates entering the proof are 
uniform in y. This implies the estimate (14. 3|) . □ 

Next, we consider scaling in the opposite direction, that is a — > 0. 
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Theorem 4.2. If u G C n+1 (IR) is periodic and \u'(y) \ + \u"(y) \ H — ■ + |n( n )(?/)| > /or some n 
and a// n, then there is C > sitc/i t/iat /or small enough a > and |v4| > Ca~ n L, the solution 
of (|4.2|) wzia initial condition (|4.1|) satisfies T(t,x,y) — > as t —+ oo, uniformly in R 2 . 

Let us give a short explanation of this result. Consider first the situation as in where 
w(n) = siny was analyzed numerically. When there is no flow, the critical quenching size, 
according to results of Kanel', is of the order £. Therefore one expects that to quench initial 
data of size L, the flow should be able to thin it down to width £, given by Theorem 13. 11 in time 
t ~ M^ 1 (before the reaction picks up). The differential of velocities near the tip at points which 
are distance ~ £ apart is Aa 2 £ 2 , so we get the condition for quenching Aa 2 £ 2 r ~ L, which is 
consistent with our theorem (since the assumptions are satisfied with n = 2 when u(y) = sinn). 
In a more general setting, assume u is smooth enough and u'(0) = u"(0) = ■ ■ ■ = n (n ~ 1 -'(0) = 
(and u' does not vanish to a higher degree elsewhere). If then A grows slower than 0(a~ n ) as 
a —* 0, the functions Au(ay) become very flat on intervals around with increasing lengths. 
The heuristic reasoning above then suggests that one should not expect quenching for small 
a's. 

To prove the theorem, we will need an auxiliary lemma. For b G iS n_1 , the unit sphere in M n , 
we define 

P b (y) = b n y n + b n ^y n - x + ■■■ + b lV . 



Lemma 4.3. Given any t > and K < oo we have 



sup f( [ P h (W Q 3 )ds e[a,a + e] 



|W°| < K for s G [0, t] ) -> 



(4.4) 



as e 



0. 



Proof. We define 

q(b,a,e) 



F 



P b {W° s )ds e[a,a + e] \W°\ < K for s G [0, t]\ 



and we let ee K n +K n ~ 1 H YK so that \P b+5 {y) - P b {y)\ < N\6\ whenever \y\<K. Hence 

we need to show that, just as p in Section [21 q — > as e — >• 0, uniformly in (6, a) G x M. 
Notice that we do not need to exclude the value tP b (0) = in the above probability because 
the P b s have no plateaux. 

The proof is identical to that of Lemma [2.21 First, the absence of plateaux in the P b s gives 
q(b, a, 0) = 0. Then with 5 = tN\8i\ + \5 2 \ + \8 3 \ we have 

\q(b + 5i,a + 5 2 ,e + 5 3 ) - q(b, a,e)\ 

< f(J^ P b (W Q s )ds G [a - 5, a + 5} U [a + e - 5, a + e + 5} \Wg \ < K for s G [0,t]^ 

which goes to zero as 5 — ► because q(b, a, 0) = q(b, a + e, 0) = 0. Thus, q is jointly continuous 
in (b,a,e). This means that q(b,a,e) — > as e — ■> 0, uniformly in any compact subset of 
(S™ -1 x E. Finally, g(6, a, e) = for |a| > tN + e, finishing the proof. □ 
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Proof of Theorem \4- 6 4 Since u G C n+1 (R) and is periodic, + \u"(y)\ + ■ — h |«( n )(?/)| > p 

for some p > and all y. Let K be such that 



PMW°| < K for s G [0,1]^ > 1 



^o e -M/2 K _ 



LetC>0, |A| > and c = Hw^H^/^+l)!. Then if b k = (ay) / A;! for /c = 1, 

Taylor's theorem gives us 

u(a(y + 5)) = u{ay) + P b {a5) + 5a n+1 |<5| n+1 

for some |c| < c. Notice that 6 need not be a unit vector here. 

With all the following probabilities conditioned by \W®\ < K for s G [0, 1], we have 



sup*( —,x,y ) 



< sup Fix — 







< 


supP 






< 


supP 




a,y 




supP 







A 
2k 



u(a{y + W° s ))ds G [-L, L] 

4ac 



a~"u(a(y + iy s u ))(is G 



a, a + 



C 



7 0_ 6 ~M/2k 



00 -M/2n 

+ 2 



P d {W° s )ds G 



4K at 

a, a + — + 
c 



+ 



^0 -M/2k 



(4.5) 



with d fc = b k a k ~ n and iV = 2ciT 1+1 . If we take a < 1, then \d\ > \b\ > p/(n + 1)! and so there 
are e G S 71 ^ 1 and r > pjin + 1)! such that d = re. The last expression in ()4.5jl is then at most 



sup P 

ee< 5n-l a 



P e {W° s )ds G 



a, a + 



[n + 1)\ /4k 



Lemma 14.31 ensures that for some C < oo and all small a the supremum is smaller than 
9 e- M/2K /2, and then flUD and (jSH) give flH3J) for r = (2/c) -1 . The result follows. □ 
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